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Abstract 

We investigate the quantum entropy, its power spectrum, and the excitation inversion of a Cooper pair box in- 
teracting with a nanomechanical resonator, the first initially prepared in its excited state, the second prepared in a 
"cat"-state. The method uses the Jaynes-Cummings model with damping, with different decay rates of the Cooper 
pair box and distinct detuning conditions, including time dependent detunings. Concerning the entropy, it is found 
that the time dependent detuning turns the entanglement more stable in comparison with previous results in literature. 
With respect to the Cooper pair box excitation inversion, while the presence of detuning destroys the its collapses and 
revivals, it is shown that with a convenient time dependent detuning one recovers such events in a nice way. 
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1. Introduction 

In the last years there has been a great interest in the production of new nonclassical states of the quantized 
electromagnetic field, one of the interesting topics of Quantum Optics. Despite the field quantization in 1925, quantum 
optical effects were observed only seven decades after, the first of them being the antibunching effect, as predicted by 
Carmichael and Walls in 1976 [ 1 ], experimentally confirmed by Kimble, M. Dagenais and L. Mandel [3 ] in 1977. A 
second nonclassical effect was observed in 1985 by Slusher et al. [4], theoretically anticipated by Stoler et al. in 
1970. A third one, the oscillations in the photon statistical distribution, was observed in 1987 by Rempe et al. (6). 
Since then, various nonclassical states of the quantized electromagnetic field were studied, including their practical 
realization in laboratories in different systems - one of them being the famous Schrodinger "cat" state, its generation 
being sugested by Yurke and Stoler (7), Davidovich et al. (8), etc; its first experimental observation was obtained by 
the group of Haroche (9). More recently, the community became aware of the first experimental observation of the 
decoherence of the Schrodinger "cat" state, in both realms of optical [10] and atomic physics iPTO . and constituting 
the first observation of the passage through the frontier that separates the quantum and classical physics. After that, 
another interesting topic emerged, as the quantum teleportation of states, first suggested by Bernnett et al. fT2l . based 
on the nonlocal character of quantum mechanics and contextualized by the EPR entangled states fT3l . This somewhat 
"bizarre" effect was first observed experimentally in 1997, by the group of Zeilinger fT4l . concerning the teleportation 
of a single photon state; the effect was later extended for atomic states and also for a huge quantity of photons fl5l . 
Then, several publications in this line appeared in the literature |[T6l[T7lfT8l[T9ll . 

Besides the nonclassical effects of light field states, many researchers became interested in the study of new states 
and new effects they could exhibit, mainly concerning with their potential applications (20). Then, it became also 
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interesting the study of various schemes for the generation of nonclassical light states (2HE2H23). To this end, two 
lines of study emerged: (i) when the issue is concerned with a state of a stationary field, inside a high-Q microwave 
cavity; (ii) when concerning with a traveling field, either thoughout the free space or a medium (optical fibers, beam 
splitters, prisms, etc.). In both cases various proposals appeared in the literature (2U|25][26][27]]. The extension of 
these investigations for atomic systems has been also implemented. In this case the system no longer concerns with 
traveling fields or a field trapped inside a high-Q cavity; instead, it consists of atoms either inside or crossing a cavity, 
including atoms inside a magneto-optical trap I28ll29l . 

When focusing either the field or the atomic case the theoretical strategy starts from a Hamiltonian describing 
the atom-field system, traditionally treated via the Jaynes-Cummings model and the atom-field coupling usually con- 
sidered as a constant parameter. Comparatively, the number of such works in the literature is very small when one 
considers the atom-field coupling and/or atomic frequency as a time dependent parameter (30j [3T] [321 |33j [34), in- 
cluding the case of time dependent amplitude [35]. Nevertheless, this scenario is also relevant; for example, the 
state of two qubits (qubits stand for quantum bits) with a desired degree of entanglement can be generated via a time 
dependent atom-field coupling l36l ; such coupling can modify the dynamical properties of the atom and the field, 
with transitions that involve a large number of photons lf3VTl . In general, these studies are simplified by neglecting the 
atomic decay from an excited level. Theoretical treatments taking into account this complication of the real world 
also employs the Jaynes-Cummings model. In these case, as expected, one finds decoherence of the state describing 
the system, since the presence of dissipation destroys the state of a system as time flows. 

Here, taking advantage of what we have learned on the atom-field interaction, we will study an advantageous 
system in practice (due to its rapid response and better controllability l38l ) by considering a nanomechanical resonator 
(NR) interacting with a Cooper pair box (CPB). This nanodevice has its own interest since its macroscopic nature and 
peculiar effects of low-frequency noise in the solid-state impose obstacles requiring more careful studies than a mere 
translation from quantum optics. It has been explored in the study of quantum nondemolition measurements 13911401 , 
in the study of decoherence of nonclassical states, as Fock states and superposition or entangled states describing 
mesoscopic systems BTlL etc. The fast advance in the tecnique of fabrication in nanotecnology implied great interest 
in the study of the NR system in view of its potential modern applications, as a sensor, largely used in various domains, 
as in biology, astronomy, quantum computation, and more recently in quantum information l42l to implement the 
quantum qubit |43 ] and in the production of nonclassical states, e.g.: Fock states l44l . Schrodinger's "cat" states l45ll , 
squeezed states |46], clusters states l47lL etc. In particular, when accompanied by superconducting charge qubits, the 
NR has been used to prepare entangled states l48l . Zhou et al. l46l have proposed a scheme to prepare squeezed states 
using a NR coupled to a CPB qubit; in this proposal the NR-CPB coupling is under an external control while the 
connection between these two interacting subsystems play an important role in quantum computation. Such a control 
is achieved via convenient change of system parameters, which can set "on" and "off' the interaction between the NR 
and the CPB, on demand. 

One of the desired goals in this report is to verify the behavior and properties of an entangled state describing 
the CPB-NR system, via the Jaynes-Cummings model, by considering the energy dissipation in the CPB during its 
transitions from an excited level to a ground state. Another target is to verify if, and in which way, the time dependence 
of the CPB-NR coupling modifies the dynamical properties of the state describing a subsystem. We will also study the 
time evolution of the quantum entropy and its power spectrum, as well as the CPB excitation inversion. There are some 
evidences of entropy production, including the fact that the power spectrum of stationary systems and subsystems can 
be used as dynamical criteria for quantum caos |49l[50]|. For the entropy power spectrum, such criteria embody those 
already discussed in the literature concerned with fixed parameters. Then, it seems adequate to look at the various 
characteristics of the entropy to formulate a reasonable and suficient universal dynamical criterium for the quantum 
caos. The degree of entanglement, represented by the entropy in certain circumstances, has also shown itself being 
sensible to the presence of a classical caos l5T][52ll . 

2. Model hamiltonian for the CPB-NR system 

There exist in the literature a large number of devices using the SQUID-base, where the CPB charge qubit consists 
of two superconducting Josephson junctions in a loop. In the present model a CPB is coupled to a NR as shown 
in Fig. ([T]); the scheme is inspired in the works by Jie-Qiao Liao et al. 1431 and Zhou et al. l46ll where we have 
substituted each Josephson junction by two of them. This creates a new configuration including a third loop. A 
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superconducting CPB charge qubit is adjusted via a voltage V\ at the system input and a capacitance C\ . We want the 
scheme ataining an efficient tunneling effect for the Josephson energy. In Fig.Q we observe three loops: one great 
loop between two small ones. This makes it easier controlling the external parameters of the system since the control 
mechanism includes the input voltage V\ plus three external fluxes ®(£), O(r) and ® e (0- In this way one can induce 
small neighboring loops. The great loop contains the NR and its effective area in the center of the apparatus changes 
as the NR oscillates, which creates an external flux ® e (0 that provides the CPB-NR coupling to the system. In this 




Figure 1: Model for the CPB-NMR coupling. 



work we will assume the four Josephson junctions being identical, with the same Josephson energy Ej, the same being 
assumed for the external fluxes O(^) and O(r), i.e., with same magnitude, but opposite sign: = -O(r) = O(x). In 
this way, we can write the Hamiltonian describing the entire system as 

H = coa^a + 4E C (n § - ^ j fr z - 4£° cos |^ j cos |^ j & X9 (1) 

where a\a) is the creation (annihilation) operator for the excitation in the NR, corresponding with the frequency oj 
and mass m\ Ej and E c are respectively the energy of each Josephson junction and the charge energy of a single 
electron; C\ and Cj stand for the input capacitance and the capacitance of each Josephson tunel, respectively. O = 
h/2e is the quantum flux and N\ = C\V\/2e is the charge number in the input with the input voltage V\. We have 
used the Pauli matrices to describe our system operators, where the states \g) and \e) (or and 1) represent the 
number of extra Cooper pairs in the superconduting island. We have: & z = \g) (g\ - \e) (e\, & x = \g) (e\ - \e) (g\ and 
E c = e 2 l(c l +ACf ) . 

The magnectic flux can be written as the sum of two terms, 

<D, = + Blx , (2) 

where the first term is the induced flux, corresponding to the equilibrium position of the NR and the second term 
describes the contribution due to the vibration of the NR; B represents the magnectic field created in the loop. We 
have assumed the displacement x described as x = xo(a^ + a), where xq = ^mco/2 is the amplitude of the oscillation. 
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Substituting the Eq.(|2| in Eq.jlj and controlling the flux we can adjust cos (^) = to obtain 

H = coa^a + 4E C \N g " ^j<^z " 4£° cos y^A sin ( j & x (3) 

and making the approximation nBixl&o « 1 we find 

H = ojcC a + -C0()cr z + Ao(a f + a)fr x , (4) 

where the constant coupling Ao = cos (^) an ^ me effective energy = 8£ c (a^ - ^ . In the rotating 

wave approximation the above Hamiltonian results as 

H = ojci^a + -cl>o& z + A()(&+a + ci f &-). (5) 

Next, we will consider a more general scenario by substituting oj — > = a) + f(t) and /io — > ^(0 = 

/lo [1 + f(t) /oj] (371 E3; in addition, we assume the presence of a constant decay rate y in the CPB, from its ex- 
cited level to the ground state; ojo is the transition frequency of the CPB and Ao stands for the CPB-NR coupling. 
<t+ and fr z are the CPB transition and excitation inversion operators, respectively; they act on the Hilbert space of 
atomic states and satisfy the commutation relations [&+,&-] = fr z and [<x z , &+] = ±&+. As well known, the coupling 
parameter A(t) is proportional to (t) IV it), where the time dependent quantization volume V (t) takes the form 
V(t) - Vol [1 + f(t) /oj] |54j[33l|53l. Accordingly, we obtain the new (non hermitean) Hamiltonian 

t 1 . y 

H = oj{t)a a + -ojod-y + A(t)(cr + a + a &-) - i- \e) (e\ . (6) 

Non hermitean Hamiltonians (NHH) have been largely used in the literature. As some few examples we mention: 
Ref. [55 1, where the authors use a NHH and an algorithm to generalize the conventional theory; Ref. 1551 , using 
a NHH to get information about entrance and exit channels; Ref. [ 57 1, using non hermitean techniques to study 
canonical transformations in quantum mechanics; Ref. 1581 , solving quantum master equations in terms of NHH; 
Ref. 1 59], using a new approach for NHH to study the spectral density of weak H-bonds involving damping; Ref. 
1551 , studing NHH with real eighenvalues; Ref. (5T1L using a canonical formulation to study dissipative mechanics 
exhibing complex eigenvalues; Ref. |62], studing NHH in non commutative space, and more recently: Ref. 1531 , 
studing the optical realization of relativistic NHH; Ref. [ 33 ], studing the evolution of entropy of atom-field interation; 
Ref. (32), using a damping JC-Model to study entanglement between two atoms, each one inside distinct cavities 



3. Solving the CPB-NR system 

Now, the state describing our time dependent system can be written as 

I* W> = S C ^n (0 IS, n) + C e , n (0 \e, n)). (7) 

Taking the CPB initially prepared in its excited state \e) and the NR in a superposition of two coherent states, \/3) = 
n(\a)+\-a)), and expanding each coherent state component in the Fock's basis, i.e., \a) = exp(-\a\ 2 /2) J]™ =0 (a n / ^/nl)\n), 
we have \j3) = ^™ =0 l^> ,where n = [2 + 2exp(-2a 2 )]" 1/2 is the normalization factor. Assuming the NR and CPB 

decoupled at t = and the initial conditions C g ^ n (0) = and Yun=o \Ce,n (0)| 2 = 1 we may write the Eq. ^ as 

l^(0)>= Y°° F n \e,n). (8) 
The time dependent Schrodinger equation for the present system is 

.d\y(t)) 



dt 



= H\*¥(t)), (9) 
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with the Hamiltonian H given in Eq. Substituting Eq.|6]) in Eq.{9} we get the (coupled) equations of motion for 
the probabilitity amplitudes C e ^(t) and C g , n+ i(t): 

= -ina>(t)C e , n {t) - l -co c e , n (t) - u(t) ^fcTTc g , n+l (t) - y -c e , n (t\ (10) 

dCa n+ \(t) i J 

— ^ = -i(n + l)a)(t)C g , n+1 (t) + -a) C g , n+1 (t) - iA(t) Sn + lC e , n (t). (11) 

The solutions of the coefficientes C e , n (f), C g , n+ \(t) furnish the quantum dynamical properties of the system, includ- 
ing the CPB-NR entanglement. 

For the cases fit) = and f(t) = const, the Eq.([TQ} and Eq.( 1J_ ) are exactly soluble. We find, analytically, 



C g , n+1 (t) = (l/Ol-lUe- 1 ' 4 * (e 1 ^ - e-W) VJ^TTfJ, (12) 
C e , n (t) = (1/20 [ie~ im (e lW iiy + 2co- if - 2oj ) - iiy + 2co + if - 2oj )) F n ], (13) 

where 6 = y + 2ioj(l + 2ri) and f = [y(y + 4i(ojQ - oj)) - 4(oj 2 + cufy - 16A 2 (1 + n) + 86ixjo] 1/2 . However, when the 
coupling f{t) is time dependent the solution of this system of equations is found only numerically. 

As well known, in the presence of decay rate y in the CPB the state of the whole CPB-NR system becomes mixed. 
In this case its description requires the use of the density operator pcv, which describes the entire system. To obtain 
the reduced density matrix describing the CPB (NR) sub-system we must trace over variables of the NR (CPB) 
sub-system. For example, 

Pnr = Tr cpb(Pcn) = YjTj + C g/l (t)CT g/ ,(f)] \n) (nr\ . (14) 

n n' 

4. Entropy of sub-systems 

Recently, researchers have employed several methods to study the dynamical of entanglement I32ll33ll34ll64ll65l . 
As proved by Phoenix and Knight [ 66 ] the von Neumann entropy offers a quantitative measure of disorder of a system 
and of the purity of a quantum state. Such entropy, defined as Snr(C) = -Triple) In Pn(C)), is a measure that is 
sensible to quantum entanglement of two interacting subsystems. The quantum dynamics described by the Eq. ^ 
furnishes the CPB-NR entanglement and we will employ the von Neumann quantum entropy as a measure of the 
degree of entanglement. The entropy S of a quantum system, when composed of two subsystems, obeys a theorem 
due to Araki and Lieb, which stablishes that: \Scpb ~ S nr\ ^ S < S qpb + S nr\ S cpb and S nr standing for the 
entropies of the subsystems. S stands for the total entropy of CPB-NR system. One immediate consequence of the 
above inequality is that, if one prepares the entire system in a pure state at t = 0, then both components of the whole 
system have the same entropy for the subsequent time evolution. So, when assuming our system initially in a pure and 
decoupled state the entropies of the CPB and NR become identical, namely, S cpbO) = SnrO)- Then, one only needs 
to calculate the quantum entropy of a subsystem to get its entanglement evolution. We obtain, from the Eqs. ( [T4] ) and 
Snr = -Trip NR lnp NR ), 

SnrO) = - [A + NR (t) ln(A^(0) + A~ NR (t) ln(A"^(0)] , (15) 

where, 



Ul± yjm\Ri) ~ {Ri\Ri)) 2 + 4|(/? 1 |7? 2 )| 2 J, (16) 



with {R X \R X ) = ZZo Knittf > <*2l*2> = Z~ =0 l^wf and (R^) = (R^T = ZZo C* e/l+1 (t)C g/l+1 (t). 

We can now look at the time evolution of the NR entropy. We will assume the NR subsystem initially in an even 
"Schrodinger-cat" state. Firstly we consider the resonant case (f(t) = 0); the time evolution of the NR entropy with 
different decay rates y in the CPB, with oj = a>o = 200(Uo and the "cat"-state with a = 5, as shown in Figs. |5Ja),|2jb), 
and[2jc). In an ideal case the CPB decay rate vanishes. As displayed in Fig. [2ja) the maximum value of the entropy of 
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the NR is close to In 2. Just after the start of the CPB-NR interaction the entropy of the NR stabilizes at this value by 
a small interval and then recovers the oscillations as time goes on. The CPB sub-system is stable while standing in its 
ground state \g); but when lying in its excited state \e) , various factors as spontaneous emission among others, imply 
its decay to the ground state. For a small decay rate (see Fig. [2fb) the maximum entanglement becomes significant 
only for large times. However, the increase of the decay rate produces a drastic change on the entanglement (see Fig. 
[2jc), with a great reduction in their swings, leading the CPB to its ground state (zero entropy). This effect upon the 
entropy of the CPB also affects the entropy of the NR (Figs. [2]). Secondly, we modify the previous case by including 
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Figure 2: Time evolution of the Entropy when the NR is initially prepared in an even "cat- state", for different values of the decay rate y: (a) 
y = O.O^o, (b) y = 0.0U 0> and (c) y = 0.05A , with a = 5, oj = oj = 2000i 0> fif) = 0. 



the presence of a detuning (f(t) = A ^ 0) to verify its influence opon our interacting system. We take the decay rate as 
y - 0.05/lo, with fit) = A = const and A <*c cjq, oj. As result the entanglement remains for long time as the value of A 
increases, as we see comparing Fig. |3ja) with Fig. |2jc); this event is accompanied by a diminution of the maximum 
entropy (see Figs[3Ja) and[3Jb)). When the detuning increases, the CPB transitions decreases (cf. Figs. [9]). 




Figure 3: Same as in Fig^ now for different values of detunings (cf 
f(t) = A = const): (a) A = l(Uo and (b) A = 20Aq , with a = 5, 
oj = ojq = 200CW , y = 0.05/i . 



Figure 4: Same as in Figs, ^and^for time-dependent detunings (cf 
f(t) = csm(ojrt) ): with (a) c = 20/to e ojr = O.lAo, and (b) c = 2(Uo e 
ojf = 0.5/i() 



Thirdly, we extend the detuning to the time dependent case, assuming f(f) = c sin(coft), where c and cot are 
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parameters of amplitude and frequency modulation of the NR with the condition ojf < c <?c ojq, oj. Comparing the Fig. 
|4ja) with Fig. |3jb) we see that the sinusoidal modulation does not favor the entanglement for long time. However, the 
frequency modulation turns the maxima of entanglement greater (see Fig. [4]). Comparing Fig. |4ja) withQb) we see 
that when the frequency ojr grows the oscillations of the entropy decrease. 



5. Power Spectrum of the Entropy 

To get a better understanding of the entropy we have considered its power spectrum (PS). It consists of a frequency- 
dependent function, being real, positive, and constructed from the following Fourier transform l54l 

PS(m)=- S NR (t)exp(imt)dt, (17) 
n Jo 

where r max = /io*max stands for the maximum interaction interval in the plot S nr(0 versus Aot. 

The entropy PS is obtained from the above equation and plotted in Figs. |2ja),[2jb) and[2jc). As expected, the 
amplitude of oscillations of this PS is reduced in the presence of growing decay rates (see Figs. [5ta),[5jb) and[5jc) 
); when we add a constant detuning ( f(t) = A ^ 0) and a decay rate y - 0.05/lo we see in Figs7[6ja) and[6jb) the 
maximum value of the PS increasing when A also increases. In the case f(t) = csin(ojft) the frequency of entropy 
PS is smoothly attenuated, with a peak around Q = 0.1, as shown in Fig. [7ja). When the frequency ojt increases the 
entropy PS is rapidly attenuated (cf. Figs. [TJa) and[7Jb)). 




0,00 0,05 0,10 0,15 0,20 



Figure 5: Time evolution of the Normalized Power Spectrum when the NR is initially prepared in "cat-state", concerning with the Entropies shown 
in the plots of Fig. |2j the plot^i) refers to the plot^i), i = a, b, c. 



6. Excitation Inversion of the CPB 

The CPB excitation inversion, IcpbQ), is an important observable of two level systems. It is defined as the differ- 
ence of the probabilities of finding this system in the excited and in the ground state; for the CPB it reads, 



oo 

WO = J] \\Ce,n(t)f ~ |C«,« + l«f 
n=0 



(18) 



The Eq. ( [T8] ) allows us to look at the time evolution of the CPB excitation inversion. First, we assume the resonant 
case {fit) = 0), for different values of the decay rate y, with a = 5 and oj = ojq - 200(Uo as in Fig. [8] Figs. 
[8ja), (b) and (c) exhibit identical collapse and revival, but with different amplitudes: the higher the decay rate, the 
lower the amplitude of oscillations of CPB excitation inversion. However, in the presence of a fixed detuning, with 
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0,00 0,05 0,10 0,15 0,20 
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Figure 6: Same as in Fig. U\ concerning with the Entropies shown in 
the plots of Fig. the plot^i) refers to the plot^i), i = a, b. 




0,00 0,05 0,10 0,15 0,20 

a 



Figure 7: Same as in Figs, \5\and \6\ c oncerning with the Entropies 
shown in the plots of Fig. |4j the plot^i) refers to the plot^i), i = a, 
b. 



fit) = A = const and A <?c oj , co, we see that CPB excitation inversion in Fig. |9ja) occurs only inside the interval 
t = Aot £ (15, 30); differently, in Fig. [9jb) this event occurs inside the interval r e (50, 75), with amplitude smaller 
than that in Fig. [9ta); this is the effect caused by a constant detuning upon the CPB excitation inversion. 





Rc)1 - 



50 100 150 200 



Figure 8: Time evolution of the CPB Excitation Inversion with the NR initially prepared in the even "cat-state", for various values of the decay 
rate: (a) y - 0.0/lo (b) y - 0.01/lo, and (c) y = 0.05/io, with a = 5, u = qjq = 2000/10, and fit) - [resonance). 



For a time dependent detuning, f(t) = csin(ojft), the frequency of the CPB excitation inversion accompanies the 
frequency cot, as shown in Figs. [10] When we compare the Fig. [H)|b) with Fig. [8jc), we see: if ojf increases, the 



interval of collapse of the excitation inversion also increases. However, looking at Figs. [T0ja),[T0jb) we see that the 
increasing of ajr as in Fig. [lOjb) introduces equally spaced collapse intervals, accompanied by revivals modulated 
by the parameter ojf . Now, we compare the CPB excitation inversion with constant detuning (fit) = A = const) 
and with a time dependent detuning (f(t) = csin(ojft))\ looking at Figs. [9] we see the plots of excitation inversion 
showing neither collapses nor revivals, with exceptions of small regions exhibinting excitation inversion (Fig[9ja)); 
in Fig[9jb) only a single such region appears. However, when considering a time dependent detuning (Fig. 10 b)), it 
nicely restitutes those collapses and revivals that appear in the resonant case (cf. Fig. [8tc)). 
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Figure 9: as in Fig. ^for diferent values of detuning (cf. f(t) ■ 

A = const): (a) A = 10/io and (b) A = 20Ao, with a = 5, a> = qjq - 
2000^0, 7 = 0.05^0- 



Figure 10: Same as in Figs, ^and^for time-dependent detunings (cf. 
fit) - csm(ajrt) ): (a) c - 2(Uo, <^ f = 0.5/lo, f^J c = 60/lo, = 20/lo> 
vWf/z a = 5, qj = qjo = 2000/lo, 7 = 0.05 /Iq. 



7. Conclusion 

We have considered a Hamiltonian model for a CPB-NR interacting system to study Entropy, its Power Spec- 
trum and the CPB Excitation Inversion. These properties characterize the entangled state that describes this coupled 
system for various values of the parameters involved. We have included dissipation and assumed the NR initially in 
a Schrodinger "cat"-state and the CPB in excited state. We have also considered the following scenarios: (/) both 
subsystems in resonance (detuning / = 0); (ii) off-resonance, with a constant detuning (/ = A ± 0), and {Hi) with a 
time dependent detuning {fit) = csin{ajn)). The results were discussed in the previous section. Concerning with the 
entropy we see that when the NR is initially in a Schrodinger "cat"-state, the entropy lasts longer than in an atom-field 
system, with the field initially in a coherent state (cf. Ref. |33|). Concerning the Excitation Inversion, an interest- 
ing result emerges: although the presence of a constant detuning destroys the collapse and revivals of the excitation 
inversion, these effects are restituted by the action of convenient time dependent detunings - even in the presence of 
damping. It is also worth emphasizing that the presence of an external force upon the NR changes the magnetic flux 
O e (cf. Fig. [TJ, which provides the control of the parameters a>{t) and X{t). 
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